Strong Correlation Effects in the Fullerene C20 
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The smallest fullerene, dodecahedral C20, is studied using a one band Hubbard model parameter- 
ized by U/t. Results are obtained using exact diagonalization of matrices with linear dimensions as 
large as 5.7 x 10^, supplemented by quantum Monte Carlo. We report the magnetic and spectral 
properties of C20 as a function of U/t and investigate electronic pair binding. Solid forms of C20 
are studied using cluster perturbation theory and evidence is found for a metal-insulator transition 
at [/ ~ 4t. We also investigate the relevance of strong correlations to the Jahn- Teller effect in C20. 

PACS numbers: 71.10.Li, 71.20.Tx, 71.30.+h 



The smallest fullerene, C20, contains 12 carbon pen- 
tagons (and no hexagons) forming a dodecahedron in a 
perfect representation of a platonic solid. Among the 
many possible isomers of C20, it is not obvious that this 
dodecahedral fullerene cage should be the most stable 
and theoretical studies [J, 111 have reached different con- 
clusions. In addition, unlike Ceo, dodecahedral C20 is not 
spontaneously formed in condensation or cluster anneal- 
ing processes Q, and its extreme curvature and strong 
reactivity led to doubts about its stability. It therefore 
created considerable excitement when Prinzbach et al. Q 
succeeded in producing the dodecahedral fullerene isomer 
of C20 in the gas-phase. Experiments have also shown ev- 
idence for solid phases of C20 @, 0] although the crystal 
structure is still debated. Density-functional studies of 
the solid forms of C20 [1, S S, @1 have suggested differ- 
ent crystal structures with the most promising candidate, 
C20 cages connected by C atoms to form a 22 atom unit 
cell 0, 0], predicted to become superconducting upon 
doping with Na.Q A simple-cubic-like phase of C2n has 
also been speculated to become superconducting [8|]. In 
their proposal for a purely electronic mechanism for su- 
perconductivity in Ceo, Chakravarty, Gelfand and Kivel- 
son [1^ have stressed the importance of structure at the 
mesoscale (Tlj. Along with the molecular solids formed 
by Ceo , solid phases of C20 would be ideal candidates for 
this picture, and a detailed understanding of these phases 
would be of great interest. Strong correlation effects are 
likely to be very important in C20, and previous stud- 
ies [ll[3,[l,0,&[i2| have treated these correlation effects 
approximatively. Here, we show that within a Hubbard 
model description, an almost exact treatment is possible 
using a large-scale numerical approach. 

Our starting point is the one-band Hubbard Hamilto- 
nian on a single C20 molecule defined as: 
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where (cia) is an electron creation (annihilation) op- 
erator, U is the on-site Coulomb interaction and riia is 
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FIG. 1: (Color online), (a) HMO levels for neutral C20 with 
Ih symmetry, (b) The Walsh diagram showing the evolution 
of the levels with the distortion, e, towards D^d symmetry, 
(c) HMO levels for the molecule with D-^d symmetry (e = 1). 



the number of electrons on site i with spin a. Due to 
the extreme curvature of C20 we expect t to be smaller 
than typical values for Ceo while U should remain close 
to that of Ceo- Consequently, we expect that a realistic 
value for U/t valid for C20 is likely larger than the value 
of U/t ^ A [ilj used for Ceo, implying that strong cor- 
relation effects play a more crucial role in the physics of 
C20. Even though Ceo and C20 share the same symme- 
try group, Ih, non- interacting C20 is a metal while Ceo 
is an insulator. This is evident from the Hiickel molec- 
ular orbitals (HMO) shown in Fig. [T^. The highest oc- 
cupied molecular orbital (HOMO), G„, is fourfold de- 
generate, containing 2 electrons for the neutral molecule. 
The lowest unoccupied molecular orbital (LUMO), Gg, 
is considerably higher in energy. Here we shall mainly 
be concerned with the neutral, 1- and 2-electron doped 
molecule which only involves the G„ levels. 

Due to the orbital degeneracy of the neutral C20 
molecule, the Jahn- Teller effect is likely to be important. 
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TABLE I: Lowest energy levels (in units of t) of the neu- 
tral C20 molecule for U — 2t,5t, labeled by spin and pseudo- 
angular momentum, for the Ih {j — jio) and Dsd [j — ja) 
configurations. 



The electronic Gu levels can couple to the Ag, Gg and 



Hg Jahn- Teller phonon modes. Theoretical studies have 
argued for the resulting lowered symmetry to be C2 [l^l , 
D5d C2h M and D^a S Here we 

follow Yamamoto et al. [l7| and assume a static defor- 
mation of D^d symmetr y. T he bond lengths for the op- 
timal Dsd structure are [l7| Uab = 1.464A, Ubc = 1.469A, 
acc' = I.519A, andocc" = 1.435A. (See Fig. 1 of Ref.lB.) 
We parameterize the distortion by letting ta, the hop- 
ping along the bond of length a, depend on a parameter 
e in the following way: ta/t = 1 — e{a — a) /a. Here, 
a — I.4712A is the average over C20 of the above bond 
lengths. Then e = and 1 correspond to the Ih and 
D^d structures, respectively. At £ = 1 the maximal de- 
viation of ta/t from 1 is less than 3.5%. In Fig. [1] the 
Walsh diagram for the evolution of the //; levels with the 
distortion e is shown along with the HMO levels for the 
optimal Dzd structure for e = 1. 

Numerics - Our exact diagonalization (ED) work is 
performed in a completely parallel fashion using SHAR- 
Cnet facilities. In addition to total particle number and 
total Sz component of the spin, we also exploit the Sio 
(Se) sub-group symmetry present in the Ih {D^d) sym- 
metry group. The basic element of this sub-group is 
a rotation of 27r/10 (27r/6) around the face of a pen- 
tagon (around a vertex) combined with reflection. We de- 
note the corresponding pseudo-angular momenta by jio 
(je). After symmetry reductions, the size of the Hilbert 
space at half-filling for the singlet states of the neutral 
molecule is 3,418,725,024 for the dodecahedral Ih con- 
figuration and 5, 699, 353, 088 for the D-^d distorted con- 
figuration. Using 64 cpu's, a Lanczos iteration for the 
Ih (Dsd) configuration is completed in 540 (980) sec- 
onds [18]. Dynamical pro perties are calculated using 
standard ED techniques [l9| . Our quantum Monte Carlo 
(QMC) work follows standard methods with ground- 
state energies obtained at T = using projector QMC 
while spectral functions are obtained from finite temper- 
ature, /3 = 10 /t, QMC (20j combined with Maximum 
Entropy methods [2l| . 

Magnetic Properties - From the non-interacting HMO 



levels in Fig. [U it would seem likely that the ground 
state of the neutral molecule is magnetic at small U/t. 
Our ED and QMC work confirms that this is the case 
for the Ih configuration for U/t < 3, where the ground 
state is observed to be an orbitally degenerate triplet, 
5 = 1, occurring at jio — 0,±2. Table |T] gives the few 
lowest energy levels, labelled by spin and pseudo-angular 
momentum, for the cases oi U/t = 2 and 5. For U/t — 5, 
we find that the ground-state for the Ih configuration is 
a non-degenerate singlet, 5 = 0, occurring at jio — 5, 
and separated from the lowest lying excitation, another 
singlet at jio = 0, ±2, ±4, by a gap of O.lt. The lowest 
triplet excitation is found at jio = 0, ±2. This ordering 
of levels continues to hold for larger U/t, although the en- 
ergy scale decreases with increasing U. The degeneracies 
and excitation gaps at large U/t agree with ED studies of 
the dodecahedral 5 = 1/2 antiferromagnetic Heisenberg 
model [i^ . In fact the ground state energy calculated for 
neutral C20 in the large U limit {U > 50) can be related 
to the ED result for the Heisenberg model to the accu- 
racy that the latter has been calculated [2^ . From Table 
H] it is clear that the system crosses over from a triplet 
to singlet ground state between U/t = 2 and U/t — 5. 
Assuming an approximately linear dependence onU/t oi 
the energy of the triplet states at jio = 0, ±2 and the sin- 
glet at jio = 5, we estimate that this transition occurs at 
Uc/t ^ 4.10, indicated by the solid vertical line in Fig. [51 
The fact that the ground-state for the neutral molecule 
for U > Uc IS a non-degenerate singlet implies that the 
molecule is stable against Jahn- Teller distortions. 

Surprisingly, the D^d distorted molecule follows the 
same pattern with the exception that at U/t = the 
unique ground-state is a singlet. However, once U/t be- 
comes of order of the splitting of the Gu levels, ^ 0.0686, 
the ground-state becomes a triplet. At U/t — 0.5 and 
2, we find that this ground-state triplet occurs at jg = 
with a number of low- lying triplet states above it. The 
Jahn- Teller distortion has therefore completely removed 
the orbital degeneracy leaving only a Kramer's degener- 
acy. At U/t = 2, the lowest-lying excitation is a triplet at 
jg = ±2, and the lowest-lying singlet, with a gap approx- 
imately twice as large, is at the same jg = ±2. As for 
the Ih configuration, we observe that the ground-state 
is a singlet at U/t = 5 occurring at jg = 3. An analysis 
similar to the Ih case yields Uc/t ^ 4.19, indicated by the 
crossed vertical line in Fig. [51 very close to the estimate 
for the Ih configuration. Again, the system remains in a 
singlet state for larger values of U /t. 

Pair Binding - The purely electronic mechanism for 



superconductivity [lOl lll | is based on a favorable pair 
binding energy. The pair binding energy is defined as 
the energy difference between having two extra electrons 
on the same and on separate molecules: 

l:\,,{N + l)^E{N + 2)~2E{N + 1) + E{N). (2) 

When negative, it is favorable to have the two electrons 
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FIG. 2: (Color online). Electronic pair binding energy 
Ab(21)/f as a function of U /t. ED (•) and QMC (o) results 
for the molecule with Ih symmetry. ED results (a) for the 
molecule with D^d (e = 1) symmetry. The solid (crossed) ver- 
tical line indicates Uc/t for neutral Ih (Dsd) molecules, (a) Aft 
versus U/t for U/t < 5. (b) At, versus t/U for 5<U/t< 100. 



on the same molecule providing a purely electronic mech- 
anism for superconductivity. Intriguingly, for neutral 
Ci2 [13 1 *s ■^^11 foi' several related models jjjl, it 
is known that is negative. For Cgo, perturbative re- 
sults indicated that the observed superconducting phase 
has its origins in a negative A;, [13, US EHl- However, 
our earlier QMC results [27| find no binding, suggesting 
that either low order perturbation theory is inadequate 
or that the QMC results are not sufficiently accurate to 
measure the binding. Thus it is of considerable interest 
to obtain exact results for C20 which can then be used to 
test the accuracy of QMC. 

Our results for Ah(21)/t are shown in Fig. [2] We first 
consider results for the Ih configuration of C2q- The 
QMC results (o) and the ED results (•) are in excel- 
lent agreement, and for U/t < 3 they show that A;, is 
positive and pair binding is suppressed. (For U > Uc it is 
not possible to perform QMC calculations due to the sign 
problem.) The ground-state for the neutral molecule is 
now a singlet and our ED results again clearly indicate 
that pair binding is not favored. As U /t increases, A^ 
reaches a minimum at U/t ~ 10. Then, as U/t — > 00, 
Afc approaches a finite positive value, consistent with the 
exact result showing the absence of pair binding for elec- 
tron doping in the U = 00 limit [28|. The perturbative 
results indicating a negative pair binding for Cgo were 
correlated with a violation of Hund's rule for two electron 
doping 10]. For C20 the same violation of Hund's rule 
occurs at large U. For U/t < 3, the ground state with 22 
electrons obeys Hund's rule and has S = 2. However, at 



U/t = 5, C20 has an intermediate spin of = 1 and at 
U/t = 6 this state has S = 0, fully violating Hund's rule. 

We can only calculate the pair binding for the Jahn- 
Teller distorted molecule in an approximate manner since 




FIG. 3: (Color online), (a) Electronic pair binding energy 
Ai,(21)/f versus the distortion, e, towards D^d symmetry, for 
for U/t = 2(0), 8(«). The inset, (b), shows the shift in the 
ground state energy of the neutral molecule, AE, versus e. 



the distortion will depend on the electron doping. We 
make the simplifying assumption that the distortion is 
static, of D^d symmetry, and independent of doping. 
Since we focus on 1 and 2 electron doping, which only 
involves Gu orbitals, it is reasonable to assume that the 
symmetry of the distortion is the same. Saito et al. [T^ 
find the same symmetry for the neutral and negatively 
charged molecule. However, the optimal e will show a de- 
pendence on the doping which we ignore. Our results for 
Afc for the D^d structure (e = 1) are shown in Fig. [2] (A). 
For all values of U/t, wc find that Af, for the D^d struc- 
ture is higher than for Ih- The Jahn- Teller distortion 
appears to work against pair binding. This is confirmed 
in Fig. [3^, where the pair binding energy is shown as a 
function of e. In Fig. we show results for the shift 
in the ground-state energy versus e where the strong de- 
pendence on e for U/t — 2 is indicative of the molecule 
being Jahn- Teller active. By contrast, for U/t — 8 the 
dependence on s is very shallow, signalling that the Ih 
structure is stable. 

Spectral Functions, Solid C20 ~ Next we calculate the 
density of states, N{ll)), and wave- vector dependent spec- 
tral functions, A{k,uj) — — (l/7r)Im[G'(fc, w + Eq + ir/)], 
for a three-dimensional solid of C20 molecules. Here, 
Eq is the ground-state energy and G the single-particle 
Green's function. The calculation is performed by clus- 
ter perturbation theory (CPT) ^ ^ using QMC and 
ED data. In all cases, delta-functions were treated as 
Lorentzians with a broadening of rj = 0.1. The QMC 
version of this method was applied earlier to the case 
of Ceo monolayers [3l|. We idealize the hypothetical 
fee C22 structure 0,13], by a model in which the bridg- 
ing C atoms are replaced by effective hopping integrals, 
t' = —t, between C20 molecules. The resulting N{uj) and 
A{k,uj), are shown in Figs. H] and [5] for U/t ~ 2 and 5. 
Also shown in Fig. [1] are the single-molecule densities of 
states. For U/t = 2, which lies in the middle of the re- 
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FIG. 4: (Color online). Density of states, N{uj), at U = 2t 
(a) from QMC and at U — 5t (b) from ED. In both panels the 
DOS for the C22 fee lattice obtained from CPT with t' = ~t 
is shown as a dashed line. 
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FIG. 5: (Color online). Spectral functions with the chemical 
potential fi determined from A'^(a;) in Fig. 3] Results for the 
C22 FCC lattice with U = 2t (a) and U = 5t (b) as obtained 
from CPT with t' = —t. The k-axis labels refer to standard 
notation for the FCC Brillouin zone. 



gion of triplet ground state, the solid of undistorted Ih 
molecules is metallic with a complicated Fermi surface, 
while, for J7/t = 5, in the region of singlet ground state, 
the solid is insulating with a gap of about lAt. 

In conclusion, we have calculated the ground state 
properties and spectral functions of the Hubbard model 
on a C20 molecule using ED and QMC. We have identified 
a ground state crossing at Uc ~ 4. It where the system 
switches from a triplet state which is unstable against a 
Jahn- Teller distortion to a gapped singlet state which is 
stable. Extending this result using CPT, we identify a 
metal-insulator transition for the bulk solid at U = Uc- 
If the symmetric, Ih, form of C20 is found to be stable, 
a possible explanation would be that it is stabilized by 
correlations resulting from U > Uc- Additional results 



will be presented elsewhere [18 
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